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ABSTRACT 

Chandra images of the elliptical galaxy NGC 4261 have revealed an anisotropy in the spatial distribution of the 
off-nuclear X-ray sources, interpreted by Zezas et al. as evidence of an association with a young stellar population. 
Our independent analysis of archive X-ray (Chandra) and optical (INT and HST) observations confirms the anisotropy 
of the X-ray sources but conducts to a different interpretation for their origin. We find that nearly half of the X-ray 
sources are associated to a globular cluster (optical counterpart) suggesting that they are accreting low-mass X-ray 
binaries. Where color index information is available, the X-ray sources are found to reside in red (metal-rich) systems. 
The luminosity function of the X-ray sources is also consistent with the one drawn from a population of Low Mass X-ray 
Binaries. We further investigate the properties of the sample of point-like sources obtained from archival optical images 
that we suggest are good globular cluster candidates and for which we find that the projected spatial distribution is 
non-homogeneous. In addition, we show that the distributions of the optical and X-ray populations are very similar, 
which leads us to conclude that the spatial anisotropy of the X-ray sources in NGC 4261 is mostly a reflection of the 
anisotropy of the globular cluster population. 

Subject headings: binaries: close — galaxies: elliptical and lenticular, cD — globular clusters: general 
— X-rays: binaries — X-rays: galaxies 
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1. INTRODUCTION 



NGC 4261 is an early- type galaxy in the Virgo 
West Cloud, at the distance of 32 Mpc (|G^vazzie^L| 
1999). It belong s to a group of galaxies {Partial 1199ft 



Nolth eTriusI 1 1 993(1 that does not show prominent sign of 
gravitational interactions. The H-band luminosity of 
NGC 4261 is 7.9 x 10 10 L corresponding to a dynam- 
ical mass of 9.5 x 10 1 M^ (using a ma ss to light ra- 
tio of about 1.2, from lZibetti et aLll2002|) . ranking this 
galaxy among the most massive galaxies in Virgo. The 
optical spectrum is consistent with a dominant po pula- 
tion of old stars (~ 11 Gyr, iGavazzi et alJ l2002|) and 
the light distr ibution is smooth with no evidence of 
sub-s tructure l(Schweizer fc Seitzerl Il992t iColbert et al.1 
HQQ2). NGC 4261 is a FRI radio galaxy (3C 270) showing 
two radio jets emanating from a nucleus hosting a super- 
massive black hole w ith a mass of about 5 x 10 8 M© 
l)Ferrarese et al. I1996D surrounded by a dusty torus an d 
a nuclear disk Uaffe et all 119961: iFerrarese et alJ ll9961. 
The X-ray nuclear power from the underlying AGN has 
been well studied by XMM and Chandra: the X-ray lumi- 
nosity is Lx ~ 5x 10 41 erg s _1 , and the X-ray spectrum is 
a power law absorbed by a column density Nh ~ 5 x 10 20 
cm~ 2 . The source appears to be embedded in diffuse hot 
gas (fcsT ~ 0.6 keV) and shows low amplitude, rapid 
variability llChiaberge et alJl200ft iSambruna et al]l200ft 
IGliozzi et al.ll2003D . 

Thanks to the high spatial resolution of Chandra, 
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IZezas et all 120031 (ZHFM03 hereafter) have discovered 
the presence of about sixty bright off-nuclear X-ray 
sources in NGC 4261, which stand out for their large- 
scale anisotropic pattern. The anisotropy extends over 
an angle of ~ 4', corresponding to a linear dimension of 
~ 40 kpc (~ 5 e ffective core radii lar ge, since r e = 4.6 kpc 
for NGC 4261: IGavazzi et al.ll2003|) . ZHFM03 interpret 
this anisotropy as evidence that the X-ray sources are as- 
sociated with young stellar population, possibly formed 
in a recent episode of star formation triggered by some 
dynamical agent along tidal tails. In spite of any optical 
evidence of either a young population (the population in 
the n uclea r regions is as old as 15 Gyrs; see lGavazzi et all 
120031 and iGoudfrooii et al.lll994[) or of recent merging 
events, ZHFM03 suggest a high mass binary origin for the 
majority of the sources, based on the significantly higher 
efficiency of the formation of High Mass X-ray Binaries 
relative to the Low Mass ones. This interpretation makes 
NGC 4261 a rather unique example, in clear contrast 
with Chandra observations of other early- type galaxies, 
where the X-ray sources are generally distributed nearly 
isotropically, trace the dominant (old) stellar popula- 
tion, and are believed to belong t o the Low Mass X- 
ray Binary population (LMXB, e.g. iVerbunt et aLlll984t 
ISarazin et aDteOOft iKim fc F abbian dl2004|) . In addition, 
a significant fraction (~ 20% -70% ) of these LMXBs is 
associated to globular clusters (GCs: ISarazin et alJl2003t 
lAngelini et al.|200lU Kundu et alJ2002UMaccarone et alJ 
l2003HMinniti et al.ll2004j) . with a preference to GCs with 
redder color index, i.e., clusters that are more met al ric h 
l|Cote et al.lli998l ISarazin et all 12553 iWest et alJEX)^ . 
In this perspective, the case of NGC 4261 stands out as 
a remarkable exception worth of further study. 

We have therefore chosen to study its case again, in 
order to deepen our understanding of the the nature of 
the X-ray sources in external galaxies. For this rea- 
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son we have used archive observations of Chandra, of 
the Isaac Newton Telescope 1 (INT), and of the Hubble 
Space Telescope (HST) (see Table 1 for a summary of 
the dataset used) to search for optical counterparts of 
the X-ray sources, and to study their relation with the 
GC population of NGC 4261. 

The outline of the paper is as follows. In §2 we describe 
the data reduction, in the X-ray and optical bands. In 
§3 we proceed on studying the properties of the optical 
point-like sources, and we cross-correlate the optical data 
with the X-ray data. We then derive the luminosity func- 
tion of the X-ray source population. In §4 we investigate 
on the azimuthal anisotropy of the X-ray sources, and on 
that of the underlying GC candidate population. In the 
context of our new findings we discuss in §5 the nature 
of the sources and their potential in tracing the history 
of assembly of NGC 4261 in the Virgo cluster. 

2. DATA REDUCTION 
2.1. X-ray data 

NG C 4261 was observed by Chandra ifWcisskop f et alJ 
12000(1 with the ACIS-S3 CCD on 2000 May 6 for a total 
exposure of 37,210 s (see Tabled for details). The data 
were collected in 1/2 sub-array mode to avoid the pile- 
up of the central AGN source. The active field of view 
is shown in Figure ^over-plotted on an optical image of 
the galaxy. The data reduction was performed with the 
package CIA0 2 v. 3.0.2 and CALDB v. 2.2. Data were 
cleaned for high background periods due to particle in- 
duced flares and the resulting exposure time is of 28.3 
ks. Individual sources were detected using a wavelet- 
detection algorithm (wavdetect in CIAO), with the same 
parameters used by ZHFM03, in three energy bands: full 
band (0.3 - 7.0 keV), soft band (0.3 - 2.0 keV) and hard 
band (2.0 — 7.0 keV). The catalog that we obtained in- 
cludes 54 sources detected above the probability thresh- 
old of 10 -6 that implies < 1 false detection due to statis- 
tical fluctuations. The positions of the X-ray sources are 
shown in FigureEl In TableElwe give the source number 
(column 1), the right ascension and declination J2000 
(column 2 and 3). Count rates (column 4) are converted 
into unabsorbed luminosities (0.2 — 10 keV) assuming a 
power-law model with photon index 1.7 and a Galactic 
line of sight colum n density of Nh= 5.8 x 10 20 cm~ 2 
{Stark et alJll992j) . ZHFM03 do not publish a catalog of 
their X-ray sources, so a comparison of the two source 
lists is not possible. However, a visual inspection of their 
Figure 1 indicate consistency. The Chandra data relative 
to the central region of thi s gal axy has also been ana lyzed 
bv lChiaberge et all l|2003[) and lGliozzi et all l|2003fl . who 
have mostly focused on the AGN emission in the X-ray, 
radio and optic al bands. 

In particular iGliozzi et al.l l)2003t) detect 5 off- nuclear 
sources: A, B, C, E coincide with our sources 38, 33, 31, 
and 14 respectively, and D, which we do not detect (our 
source 19 is the closest, at ~ 10"). Although these au- 
thors comment of the fact that they do not find optical 
counterparts for any of their sources, they do not com- 
ment on the possible connection of these sources with the 
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jet, in spite of the evidence that two of them, (A and E) 
are aligned and at opposite ends of it. 

2.2. Optical data 

We have used archive observations of INT and HST 
WFPC2 for a systematic search of the optical counter- 
parts of the X-ray sources. 

We extracted a deep r' band image from the INT- Wide 
Field Camera (WFC) archive, consisting of three expo- 
sures of 500 seconds each, for a total integration time of 
1500 seconds. The image covers the whole galaxy (the 
field of view of the INT WFC is of about 0.5 deg 2 ) and 
has a mean seeing of - 1" l|Khosroshahi et all2 004'). The 
reduction of the images was carried out using standard 
tasks in the IRAF package 3 . Bias subtraction and flat- 
field normalization were applied to each exposure using 
the median of several bias frames and flat-field exposures 
obtained during the observations. The three exposures 
were combined using a median filter, that automatically 
removes cosmic rays. Using reference catalogues of astro- 
metric stars, we obtained a WCS solution for the whole 
field having a mean RMS of —1 arcsec. 

HST WFPC2 observations cover the central portion 
of the galaxy, as shown by the white boxes in Figure 
HI The south-east (SE) field was observed using three 
different filters: F675W, F547M, F791W while only the 
F702W filter is available for the north-west (NW) field. 
The exposure times of the four images are given in Table 
□ 

For our analysis we started from pipelined processed 
data. For both HST WFPC2 and INT, we first fit- 
ted ellipses to the isophote of the galaxy using the task 
ELLIPSE from the IRAF package, and subtracted the 
model from the data. The resulting i mage was then used 
to de tect sources with SExtractor ijBertin fe Arnoutsl 
1996). The last isophote has a major and minor axis of 
4.5 and 3.2 arcmin respectively. We improve the relative 
HST astrometry from — 0.7 to — 0.2 arcsec comparing 
the overlap region between the NW and SE field. On the 
contrary we did not try to make any improvement in the 
relative astrometry between Chandra and INT /HST im- 
ages. We note that, for both HST and INT images, the 
model subtraction leaves a quadrupole structure in the 
central region. Although this region is smaller in the HST 
data, we exclude a wider area of ~0.5 square arcmin, due 
to the presence of residuals in the INT data that could 
induce spurious results in our comparisons. Th is struc- 
ture was already noticed bv lColbert et alJ l)2001|) for this 
object, and was attributed to the boxyness of the galaxy. 

Using SExtractor we extracted point-like objects from 
the model subtracted images of both INT and HST 
datasets. We set a SExtractor detection threshold of 
2.5 a over the local background. We select as point-like 
sources objects with the following properties: (i) a mag- 
nitude error smaller than 0.1, (ii) ellipticity smaller than 
0.5, (iii) FWHM between 1 and 4 pixels (0.1-0.4 arc- 
sec) for HST, as suggested bv lGebhardt fc Kissler-Patiel 

3 IRAF is the Image Analysis and Reduction Facility made avail- 
able to the astronomical community by the National Optical As- 
tronomy Observatories, which are operated by AURA, Inc., under 
contract with the U.S. National Science Foundation. STSDAS is 
distributed by the Space Telescope Science Institute, which is oper- 
ated by the Association of Universities for Research in Astronomy 
(AURA), Inc., under NASA contract NAS 5-26555. 
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( 1999), and FWHM between 3 and 10 pixels (1-3 arcsec) 
for INT. Note that for the HST-SE field we based our 
selection upon the F675W filter owing to its longer ex- 
posure. As a final step, we have checked the SExtractor 
results by visual inspection, and we have removed objects 
on or close to the CCDs edges for the HST images. 

3. RESULTS 

3.1. Selection of candidate globular clusters 

For the SE field of the HST observation, we have been 
able to estimate the color of the sources presented in 
§2.2. The instrumental magnitudes, provided by SEx- 
tractor, were converted in the Johnson-Cousin V, R and 
I systems using the synthetic transformations given in 
Holt zman et "all 1 19951 In order to estimate the photo- 
metric error due to model subtraction we measured the 
photometry on the unsubtracted image using local back- 
ground and apertures from the subtracted image. The 
difference in the V-I color of each globular cluster ob- 
tained from the two different estimates is shown in Fig- 
ure [3J The distribution appears fairly gaussian with a 
mean difference of —0.003 ±0.07 mag (see Figure |3J). 
The total photometric uncertainty is thus given by the 
combination of the error on the measure of GCs mag- 
nitudes (given by SExtractor), the err or on the zero- 
point (taken from iHoltzman et alJll995|) and the uncer- 
tainty due to the model subtraction (0.07 mag, FigEJ. 
This results in a mean error on the V-I color of ~ 0.11 
mag. We used F547M and F791W filters to determine 
the color index (V-I) of the extracted objects. We se- 
lect 325 objects detected in e ach of the three filters with 
V-I c olor between and 2 ijGebhardt fc Kissler-Patid 
1999), that we consider GC. Their color index and mag- 
nitude distributions shown in Figure 0] are consistent 
with those of the globular cluster populations observed in 
other massive early-type galaxies (!A^h inaii^L^epJ[199^; 
[ Cote et al.ll998HKundu fc Whitmorei200lULarsen et all 
l2TMTIWest^t^lT]20"o^r " reinforcing our definition. The 
color index distribution spreads from V-I ~ 0.5 to ~ 2, 
indicating that both populations of GC (red and blue) 
are present. Thus the unimodal distribution observed in 
the V-I color histogram is probably due to the blending 
of different GCs populations. 

Since the color information is restricted to the HST- 
SE field, we cannot classify as GCs the sources extracted 
from the other optical images at the same confidence. 
However, we have first considered sources in the region 
common to all three datascts (~ 1 arcmin 2 ). The cross- 
correlation between the three lists indicates that ~ 97% 
(37/38) and - 83% (49/59) of the sources in the HST- 
NW and INT field respectively coincide with a GC in the 
HST-SE field. Given this high fraction, we are confident 
that even the point-like sources from the INT dataset are 
good "GC candidates", in spite of the lower resolution 
of the instrument. In what follows, we will use the "GC 
candidate" list (from the INT observation), that provides 
an homogeneous set of optical sources over the whole 
galaxy. 

As a final check, we have also compared the sky pro- 
jected distribution of the INT and HST GC candidates. 
Since HST-SE and HST-NW fields have different expo- 
sure times, we restrict the HST sample to a magnitude 
R = 24.5 (i.e. the magnitude limit of the shallowest 
field). The comparison is done through an adaptive ker- 



nel density analysis. This technique is a two-step proce- 
dure which first applies a pilot smoothing to estimate the 
local density of GCs, starting from the list of their po- 
sitions only, and then uses a smoothing window variable 
in size that decreases with increasing local density. In 
this way the statistical noise in low-density regions can 
be reduced witho ut smearing out the high density peaks 
l)Silvermanlll986|) . As shown in Figure [SJ the iso-density 
contours from the adaptive kernel analysis are in very 
good agreement, confirming the reliability of the INT 
GCs dataset. What is apparent in the figure is the pe- 
culiar distribution of the GC population, that shows two 
main concentrations NE and SW of the nuclear region. 
The presence of the same structures in both datasets en- 
sures that it is not a feature derived from the poor sub- 
traction in the central region (the black area in FigureEJ 
since the equivalent area in the HST data is negligible. 

3.2. Cross-correlation of X-Ray and Optical data 

In order to investigate the nature of the X- 
ray sources we have looked for optical counterparts. 
The c entral source (n . 29 ) is identified with the 
AON llChiaberee et all 120031 ISambruna et~aTl 120031 
iGliozzi et al.ll2003j) . and sources no. 27, 43 and 53 with 
background/foreground galaxies. The positions of the 
remaining X-ray sources have been cross-correlated with 
the GC sample, from the HST fields, with a matching 
radius of 0.5 arcsec, and from the INT dataset, with a 
matching radius of 1 arcsec. These choices take into ac- 
count only the intrinsic astrometric uncertainties of the 
three datasets (see sections 2.1 and 2.2) and not the rela- 
tive uncertainties eventually present between the optical 
and X-ray observations. Source 54 falls outside the model 
subtracted area, while sources no. 19, 31 and 33 fall in 
the central region excluded from our previous analysis. 
We found an optical counterpart for 23 of the 46 discrete 
X-ray sources that we assume all belonging to NGC 4261. 
The identification is indicated in Tabled To calculate 
the chance coincidence probability we performed 1000 
Monte-Carlo simulations by shifting the positions of the 
INT and HST optical samples by a random offset and 
performing again the cross-correlation. We adopted a 
shift larger than the search radius but small enough so 
the sources do not fall outside the distribution of the X- 
ray sources (~1.5 arcmin). The probability of finding 
more than one spurious association is ~2.5% and ~0.1% 
in the INT and HST sample, respectively. In Figure H3 
we give a few examples that show the positions of the 
X-ray sources over-plotted onto the HST SE, HST NW 
and INT fields. The figure also demonstrates the good 
correspondence between the HST and INT counterparts, 
indicating again that the identifications obtained in the 
INT field not covered by HST are not severely hampered 
by the lower resolution of the INT image. Note that, 
where the color information is available, we found that 
all the X-ray sources inhabit a red (color index V-I > 1) 
GC, see Figure 0J The optical identification with GC 
candidates for ~ 50% of the X-ray sources su ggests a 
LMXB origin for them fVerbunt fc Lewinll2004D . 

3.3. Luminosity Function of X-ray sources 

We verified the above statement by comparing the cu- 
mulative X-ray luminosity function (XLF) of our sources 
with those derived from observations of other galaxies 
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(Figure [71) . In particular, we use the funct i onal forms 
derived bv lGrimm etaD (|2003l) and lGilfanovl (|2004fl con- 
sidered by these authors as a signature of a universal 
behavior for the HMXB and the LMXB populations, re- 
spectively; namely a single power-law for HMXB, with 
slope of a = —0.61, and a broken power-law for LMXB, 
with slopes of a.\ = —0.8, and a-i = —4 for the high 
luminosity end 4 . In testing the consistency with these 
universal functions, we have derived the XLF using all 
sources in Tabled except for the AGN (source no. 29) 
and the 3 interlopers (sources no. 27, 43 and 53). We 
performed a log-log maximum likelihood fitting proce- 
dure for both HMXB and LMXB functional forms as- 
suming fixed slopes. We considered as free parameter 
the normalization factors and, in the case of a broken 
power law, the luminosity of the break. The resulting 
XLFs are shown in Figure [7J It is unquestionable that a 
single power-law cannot reproduce the LF of the sources 
in NGC 4261. It is also apparent that the LMXB broken- 
power law is a good representation of the data given also 
the value of the break luminosity, that we have left as a 
free parameter in the fitting procedure, and that is con- 
sistent with the Gilfanov value of Lb = 5.1 • 10 38 ergs~ 1 . 
This result is consistent with the hypothesis that the X- 
ray sources in this galaxy are associated to a population 
of LMXBs. Note that we do not take into account incom - 
pleteness at low luminosities ijKim fc Fabbiancl 12004); 
however, we point out that at least at the high luminosity 
end (where incompleteness should not significantly affect 
the XLF slope) the luminosity function is much steeper 
than the XLF of the HMXBs. 

4. AZIMUTHAL ANISOTROPY IN THE GLOBULAR 
CLUSTER DISTRIBUTION 

ZHFM03 have shown with a Kolmogorov-Smirnov test 
that the spatial distribution of the X-ray sources in 
NGC 4261 is anisotropic. Using our list of 50 sources 
(sources n. 29, 27, 43 and 53 excluded) we confirm the 
spatial anisotropy over a scale of about 4' with very high 
confidence (P> 99.9% from the Kolmogorov-Smirnov 
test). In the preceding sections we have shown that the 
properties of the X-ray sources in NGC 4261 are consis- 
tent with a population of LMXBs. This however does 
not explain yet their spatial anisotropy, which is in stark 
contrast with the smooth isophotal distribution of the op- 
tical light. Given the close association with the GC pop- 
ulation, which we find also has a peculiar distribution, 
we have considered the possibility that the anisotropy 
is nothing but a manifestation of an anisotropy seeded 
in the GC population. In order to test this hypothesis 
we have first done a statistical analysis on the spatial 
distribution of the GC candidates. 

In order to investigate both the details of the spatial 
distributions of the X-ray and optical sources separately 
and their relation we have compared the normalized sur- 
face density of sources again with an adaptive kernel den- 
sity analysis on the two sets of sources. In Figure [8] we 
show the iso-intensity contours extracted from the kernel 
analysis on the GC candidate sample (thin line) and on 
the X-ray sources (thick line). With this figure as a guide, 
we have defined 6 regions (see Figure [5J right panel) in 
which we have simply computed the surface number den- 

4 These values are taken from Table 3 of lGilfanovl Oool) . 



sity of sources for both GCs and X-ray sources, which 
plot in Figure El (left panel) as a function of the region 
number. The errors associated with each entry reflects 
the total number of sources in the region, and gives a 
direct measure of the significance of the observed region- 
to-region variations, We also test the correlation between 
the two datasets using the IRAF task CR0SSC0R and we 
find that the highest correlation between the two datasets 
is consistent with a null shift. 

5. DISCUSSION AND CONCLUSIONS 

Our combined analysis of the X-ray and optical ma- 
terial available for NGC 4261 shows that the large ma- 
jority of the X-ray sources in NGC 4261 are accreting 
LMXBs. About 50% of these sources have an optical 
counterpart, that we identify with a GC. Where color in- 
formation is available, the sources are associated with the 
red (metal-rich) population, as already observed in sev- 
eral other galaxies. Further evidence of a LMXB origin 
comes from the interpretation of the luminosity function, 
consistent with the "universal" distribution proposed by 
Gilfan^(|2003). 

The photometric analysis of the GC population of 
NGC 4261 carried on a portion of the galaxy (HST-SE 
field) shows that the properties of the GCs ar e similar to 
those observed in other early-type galaxies llCote et alJ 
[l998t IGebhardt fc Kissler-Pati Jl99a IWest et al.1 12004) . 
Furthermore this similarity extends to the connection be- 
tween X-ray sources and GC population ijSarazin et alJ 
2003). In particular we note that the X-ray sources are 
preferentially associated with the brightest and reddest 
GCs. As in other ellipticals, we find that the GCs host- 
ing a X-ray source are ~7% and that the X-ray sources 
associ ated with GCs are ^50%, consistent with the re- 
sult of iKundu et al J [20021 It is therefore not surprising 
that NGC 4261 ranks among those galaxies having a siz- 
able fraction of X-ray sources in GCs. The peculiarity of 
NGC 4261 resides mainly in the large scale anisotropy of 
the X-ray sources. Our adaptive kernel density analysis 
has highlighted the presence of a peculiar spatial distri- 
bution in the GC candidate population, and the existence 
of a correlation between the spatial anisotropy of the X- 
ray sources and the main over-densities seeded in the GC 
system. To our knowledge, this is the first time that X- 
ray sources indicate the presence of substructures in the 
distribution of GC. 

A non uniform spatial distribution of the GC system 
of NGC 4261 may reflect a peculiar history of formation 
of this galaxy. In the paradigm of hierarchical cluster- 
ing cosmogony, elliptical galaxies may be the end result 
of binary or multiple major mergers among gas-rich pro- 
genitors. Major mergers likely redistribute the original 
GC systems of the progenitor galaxies and trigger the 
formation of new GCs in the center as well as along the 
tidal tails of the new object arising from the merger. 
This would result in the coexistence of two populations, 
metal rich and metal poor clusters, that have different 
photometric properties, spatial distribution and dynam- 
ical ch aracteristics, as shown in a recent work bv lLi et alJ 
( 2004) . Hierarchical build-up also implies a second way 
by which galaxy can grow, namely accretion of dwarf 
satellites that would lose their GCs to the primary sys- 
tem. GCs coming from shredded satellites could also 
contribute to spatial anisotropies and would be mostly 
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metal-poor if such satellites are similar to nearby dwarfs. 
Evidence that either of these two mechanisms can ac- 
count for a lar ge fraction of the GC systems in ellipticals 
is still lacking l)Strader e t al. 2004|). A close inspection of 
the photometric and kinematical properties of the GCs 
can in principle help disentangle many aspects of such 
complex galaxy genesis. NGC 4261 could be an excellent 
target for further investigations: combined Chandra and 
HST observations can thus potentially shed light on how 
early-type galaxies were assembled. In parallel, we are 
undertaking a theoretical study running high-resolution 
simulations of galaxy mergers similar to those carried on 
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Fig. 1 — R band image of NGC 4261 from the INT telescope. The active field of the Chandra ACIS-S3 CCD and HST WFPC2 fields 
are superposed in black and white, respectively. 
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Fig. 2. — Chandra S3 image of NGC 4261. Open squares correspond to X-ray sources with optical counterpart while open circles to 
sources without optical counterpart. Excluded from the analysis are the X-ray sources (diamonds) that fall in the region where the model 
subtraction is too noisy, and source no. 54 located outside the galaxy radius. The black cross corresponds to the AGN. "Plus" symbols 
correspond to background/foreground galaxies. 
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Fig. 3. — Histogram showing the difference in the estimate of the V-I color between unsubtracted and model subtracted HST images. 
The distribution of the GCs hosting an X-ray source is indicated in black. 
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Fig. 4. — Left panel: Color index histogram of GCs in HST SE field. Right panel: Magnitude histogram for the GCs in HST SE field. In 
both panels the distributions of the GCs hosting an X-ray source are superposed in grey. Note that the X-ray sources inhabit preferentially 
metal rich (red), brighter GCs. 
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Fig. 5. — Adaptive kernel density analysis: Iso-density contours (thin line) of the analysis carried on the sample of GC candidates from 
INT are superposed to the iso-density contours (thick line) of the GC sample for HST fields. The first contour represents the 2 a level. 
Black and white regions are excluded from the extraction of the optical sources. The two main overdensities coincide spatially in both data 
sets. 



11 



FlG. 6. — Selected X-ray sources listed in Tableland indicated with a star in column 1. The underlying images are from the HST SE 
(left), HST NW (middle) and INT (right) fields, respectively. Black circles (of radius equal to 0.7 arcsec for HST and 1 arcsec for INT) 
indicate the combined astrometric uncertainties in the position of each Chandra source. 
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Fig. 7. — The integral luminosity function in the 0.2-10 keV band of the X-ray sources in NGC 4261. The solid line gives a broken 
power-law with fixed slopes, from Gilfanov (2004), rcnormalizcd to our data. We overlay the single power-law (dashed line) with a slope 
o?3 = —0.61 as suggested by Grimm et al. (2003) to describe a population of HMXBs. 
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Fig. 8. — Iso-density contours from our adaptive kernel density analysis on the GC candidates from INT (thin line) and on the X-ray 
sources (thick line). The first contour represents the 2 a level. Note that the main over-densities in the distribution of X-ray sources 
coincide with those of GC candidates. 
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Fig. 9. — Left panel: Histogram of the surface number density (number of objects per square arcmin) of GC candidate s (solid li ne) and 
of X-ray sources (dashed) as a function of their position in the field. Error bars are ±lcr computed according to IGehrelsl fl986). Right 
panel: The regions corresponding to each bin are indicated by squares superposed to the INT-model subtracted image. 



TABLE 1 

Journal of all observations used in this work. SMAG indicates the correspondent Johnson-Cousin filter. 



Observation 


P.I. 


date 


instrument 


filters 


SMAG 


exposure time (s) 


optical 


Ford 


13 Dec 1995 


HST 


F675W 


R 


2000 


optical 


Ford 


13 Dec 1995 


HST 


F547M 


V 


800 


optical 


Ford 


13 Dec 1995 


HST 


F791W 


I 


800 


optical 


Sparks 


04 May 1995 


HST 


F702W 


R 


280 


optical 


Osmond 


04 - 10 Feb 2000 


INT 


r' 


R 


1500 


X-ray 


Birkinshaw 


06 May 2000 


Chandra 


none 




37210 
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TABLE 2 

List of sources detected in the Field of NG C4261. Count rates are in the full 0.3-7.0 keV band as given by wavdetect. 

X-RAY LUMINOSITIES ARE CALCULATED ASSUMING A POWER-LAW SPECTRUM WITH AN INDEX OF 1.7 AND A DISTANCE OF 32 MPC. CPT 
INDICATES THE ASSOCIATION WITH A GLOBULAR CLUSTER (GC) IN THE HST FIELD, A CANDIDATE GLOBULAR CLUSTER (GCC) FROM THE 
INT LIST, OR THE NAME OF THE FOREGROUND/BACKGROUND GALAXY. THE QUESTION MARK INDICATES THAT NO OPTICAL COVERAGE IS 

AVAILABLE. 



source 
number 




R.A. 

(h m s) 




Dec 
(h m s) 


CountRate L x (0.2 - WkeV ) 
(xlO- 4 ^ 1 ) (xlO'^s- 1 ) 


CPT 


1 


12 


: 19 : 


14.44 


+05 


48 


48.12 


3.76 ± 1.22 


4.93 




2* 


12 


: 19 : 


16.81 


+05 


49 


50.43 


10.65 ± 1.97 


13.96 


GCC 


3 


12 


: 19 : 


18.04 


+05 


49 


54.14 


2.68 ± 0.10 


3.52 


GCC 


4 


12 


: 19 : 


18.23 


+05 


49 


12.10 


2.03 ± 0.89 


2.67 




5 


12 


: 19 : 


18.67 


+05 


47 


43.23 


4.15 ± 1.22 


5.45 




6 


12 


: 19 : 


19.18 


+05 


48 


49.31 


2.63 ± 0.10 


3.45 


GCC 


7 


12 


: 19 : 


19.33 


+05 


49 


06.53 


2.02 ± 0.87 


2.65 




8 


12 


: 19 : 


19.78 


+05 


48 


22.33 


4.03 ± 1.22 


5.28 




9 


12 


: 19 : 


19.91 


+05 


48 


50.90 


4.21 ± 1.22 


5.52 




10 


12 


: 19 : 


20.27 


+05 


49 


08.58 


7.88 ± 1.69 


10.33 


GCC 


11 


12 


: 19 : 


20.44 


+05 


48 


02.72 


5.18 ± 1.37 


6.79 




12 


12 


: 19 : 


20.60 


+05 


48 


35.89 


4.97 ± 1.37 


6.52 


GC 


13 


12 


: 19 : 


20.66 


+05 


49 


29.10 


2.73 ± 0.10 


3.58 


GCC 


14 


12 


: 19 : 


20.93 


+05 


49 


27.09 


5.83 ± 1.46 


7.64 


GCC 


15 


12 


: 19 : 


21.00 


+05 


48 


44.07 


3.98 ± 1.22 


5.22 


GC 


16 


12 


: 19 : 


21.07 


+05 


50 


17.38 


3.26 ± 1.12 


4.27 




17 


12 


: 19 : 


21.08 


+05 


48 


07.21 


2.36 ± 0.94 


3.09 




18 


12 


: 19 : 


21.20 


+05 


47 


52.14 


5.55 ± 1.46 


7.27 


GCC 


19 


12 


: 19 : 


22.00 


+05 


49 


47.02 


3.43 ± 1.12 


4.50 


? 


20 


12 


: 19 : 


22.16 


+05 


50 


59.57 


3.07 ± 1.06 


4.02 




21 


12 


: 19 : 


22.35 


+05 


48 


47.95 


4.27 ± 1.27 


5.60 


GC 


22 


12 


: 19 : 


22.54 


+05 


50 


17.00 


5.75 ± 1.46 


7.54 


GCC 


23 


12 


: 19 : 


22.67 


+05 


47 


23.04 


3.69 ± 1.17 


4.84 




24 


12 


: 19 : 


22.81 


+05 


46 


44.62 


3.45 ± 1.12 


4.52 




25 


12 


: 19 : 


22.99 


+05 


51 


11.50 


1.60 ± 0.79 


2.10 




26 


12 


: 19 : 


23.00 


+05 


50 


58.55 


5.53 ± 1.41 


7.25 




27 


12 


: 19 : 


23.10 


+05 


47 


41.43 


7.16 ± 1.66 


9.38 


VCC 344 


28 


12 


: 19 : 


23.20 


+05 


48 


07.16 


2.69 ± 0.10 


3.53 




29 


12 


: 19 : 


23.21 


+05 


49 


29.76 


1249.41 ± 21.02 


1638.28 


AGN 


30 


12 


: 19 : 


23.23 


+05 


50 


13.44 


3.44 ± 1.12 


4.51 


GCC 


31 


12 


: 19 : 


23.50 


+05 


49 


36.06 


11.02 ± 4.07 


14.45 


? 


32 


12 


: 19 : 


23.61 


+05 


47 


55.54 


4.70 ± 1.32 


6.16 


GC 


33 


12 


: 19 : 


24.01 


+05 


49 


26.42 


9.92 ± 3.02 


13.00 


? 


34 


12 


: 19 : 


24.25 


+05 


47 


54.75 


5.92 ± 1.50 


7.76 




35 


12 


: 19 : 


24.60 


+05 


51 


04.86 


8.18 ± 1.73 


10.73 


GCC 


36 


12 


: 19 : 


24.80 


+05 


50 


57.37 


2.40 ± 0.94 


3.15 




37* 


12 


: 19 : 


24.83 


+05 


50 


08.61 


3.69 ± 1.17 


4.84 


GC 


38 


12 


: 19 : 


24.93 


+05 


49 


32.11 


5.50 ± 1.41 


7.22 


GC 


39 


12 


: 19 : 


25.50 


+05 


47 


41.85 


6.53 ± 1.54 


8.57 


GCC 


40 


12 


: 19 : 


25.55 


+05 


49 


46.17 


4.92 ± 1.41 


6.45 


GC 


41* 


12 


: 19 : 


25.65 


+05 


50 


28.75 


2.73 ± 1.06 


3.58 


GC 


42 


12 


: 19 : 


25.95 


+05 


52 


20.50 


3.96 ± 1.22 


5.19 




43* 


12 


: 19 : 


26.06 


+05 


50 


13.69 


2.76 ± 0.10 


3.62 


galaxy 


44 


12 


: 19 : 


26.89 


+05 


50 


44.15 


5.74 ± 1.46 


7.52 


GC 


45 


12 


: 19 : 


29.22 


+05 


50 


00.59 


2.91 ± 1.06 


3.82 


GCC 


46 


12 


: 19 : 


29.78 


+05 


51 


14.35 


7.46 ± 1.66 


9.78 




47 


12 


: 19 : 


29.97 


+05 


50 


28.20 


5.64 ± 1.46 


7.39 




48 


12 


: 19 : 


30.51 


+05 


50 


07.97 


3.11 ± 1.06 


4.08 




49 


12 


: 19 : 


31.08 


+05 


50 


07.96 


2.50 ± 0.10 


3.28 


GCC 


50 


12 


: 19 : 


31.98 


+05 


51 


54.46 


5.60 ± 1.46 


7.35 


GCC 


51 


12 


: 19 : 


32.08 


+05 


51 


49.47 


4.85 ± 1.41 


6.35 




52 


12 


: 19 : 


33.93 


+05 


50 


35.75 


5.65 ± 1.46 


7.41 




53 


12 


: 19 : 


35.52 


+05 


50 


48.93 


4.04 ± 1.27 


5.30 


NGC 4264 


54 


12 


: 19 : 


36.87 


+05 


51 


56.22 


4.44 ± 1.32 


5.82 
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